ARL-RP-0544  •  Sep  2015 


arl 

US  Army  Research  Laboratory 


On  Finsler  Geometry  and  Applications  in 
Mechanics:  Review  and  New  Perspectives 


by  JD  Clayton 


A  reprint  from  Advances  in  Mathematical  Physics.  Volume  2015,  Article  ID  828475,  p.  11. 
[accessed  2015  Sep  8].  http://dx.doi.org/10.1155/2015/828475. 


Approved  for  public  release;  distribution  is  unlimited. 


Report  Documentation  Page 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 


1.  REPORT  DATE 

SEP  2015 


2.  REPORT  TYPE 


4.  TITLE  AND  SUBTITLE 

On  Finsler  Geometry  and  Applications  in  Mechanics:  Review  and  New 
Perspectives 

6.  AUTHOR(S) 


7.  PEREORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

US  Army  Research  Laboratory,, ATTN:  RDRL-WMP-C„Aherdeen 
Proving  Ground„MD,  21005 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED 

21-11-2014  to  18-01-2015 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

A  reprint  from  Advances  in  Mathematical  Physics.  Volume  2015,  Article  ID  828475,  p.  11.  [accessed  2015 
Sep  8].http://dx.doi.org/10.1155/2015/828475 

14.  ABSTRACT 

In  Finsler  geometry,  each  point  of  a  base  manifold  can  be  endowed  with  coordinates  describing  its  position 
as  well  as  a  set  of  one  or  more  vectors  describing  directions,  for  example.  The  associated  metric  tensor  may 
generally  depend  on  direction  as  well  as  position,  and  a  number  of  connections  emerge  associated  with 
various  covariant  derivatives  involving  afilne  and  nonlinear  coefficients.  Finsler  geometry  encompasses 
Riemannian,  Fuclidean,  and  Minkowskian  geometries  as  special  cases  and  thus  it  affords  great  generality 
for  describing  a  number  of  phenomena  in  physics.  Here,  descriptions  of  finite  deformation  of  continuous 
media  are  of  primary  focus.  After  a  review  of  necessary  mathematical  definitions  and  derivations,  prior 
work  involving  application  of  Finsler  geometry  in  continuum  mechanics  of  solids  is  reviewed.  A  new 
theoretical  description  of  continua  with  microstructure  is  then  outlined,  merging  concepts  from  Finsler 
geometry  and  phase  field  theories  of  materials  science. 

15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

19a.  NAME  OE 

ABSTRACT 

OF  PAGES 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

17 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


NOTICES 

Disclaimers 

The  findings  in  this  report  are  not  to  be  eonstrued  as  an  offieial  Department  of  the 
Army  position  unless  so  designated  by  other  authorized  doeuments. 

Citation  of  manufaeturer’s  or  trade  names  does  not  eonstitute  an  offieial 
endorsement  or  approval  of  the  use  thereof 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


ARL-RP-0544  •  Sep  2015 


ARL 

US  Army  Research  Laboratory 


On  Finsler  Geometry  and  Applications  in 
Mechanics:  Review  and  New  Perspectives 


by  JD  Clayton 

Weapons  and  Materials  Research  Directorate,  ARL 


A  reprint  from  Advances  in  Mathematical  Physics.  Volume  2015,  Article  ID  828475,  p.  11. 
[accessed  2015  Sep  8].  http://dx.doi.org/10.1155/2015/828475. 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  the  collection  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  the 
burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid 
0MB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 


5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 


11.  SPONSOR/MONITOR'S  REPORT  NUMBER(S) 

12.  DISTRIBUTION/ AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 

A  reprint  from  Advances  in  Mathematical  Physics.  Volume  2015,  Article  ID  828475,  p.  1 1.  [accessed  2015  Sep  8]. 
http://dx.doi.org/10.1155/2015/828475. _ 

14.  ABSTRACT 

In  Finsler  geometry,  each  point  of  a  base  manifold  can  be  endowed  with  coordinates  describing  its  position  as  well  as  a  set  of 
one  or  more  vectors  describing  directions,  for  example.  The  associated  metric  tensor  may  generally  depend  on  direction  as 
well  as  position,  and  a  number  of  connections  emerge  associated  with  various  covariant  derivatives  involving  affine  and 
nonlinear  coefficients.  Finsler  geometry  encompasses  Riemannian,  Euclidean,  and  Minkowskian  geometries  as  special  cases, 
and  thus  it  affords  great  generality  for  describing  a  number  of  phenomena  in  physics.  Here,  descriptions  of  finite  deformation 
of  continuous  media  are  of  primary  focus.  After  a  review  of  necessary  mathematical  definitions  and  derivations,  prior  work 
involving  application  of  Finsler  geometry  in  continuum  mechanics  of  solids  is  reviewed.  A  new  theoretical  description  of 
continua  with  micro  structure  is  then  outlined,  merging  concepts  from  Finsler  geometry  and  phase  field  theories  of  materials 
science. 

15.  SUBJECT  TERMS 

differential  geometry,  Finsler  spaces,  continuum  mechanics,  microstructure,  phase  field,  crystals 

17.  LIMITATION  18.  NUMBER 
OF  OF 

ABSTRACT  PAGES 

UU  16 

Standard  Form  298  (Rev.  8/98) 
Prescribed  by  ANSI  Std.  Z39.18 


19a.  NAME  OF  RESPONSIBLE  PERSON 

JD  Clayton 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

410-278-6146 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

Unclassified 

Unclassified 

Unclassified 

8.  PERFORMING  ORGANIZATION  REPORT  NUMBER 

ARL-RP-0544 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 


5d.  PROJECT  NUMBER 


3.  DATES  COVERED  (From  -  To) 

21  November  2014-18  January  2015 

Sa.  CONTRACT  NUMBER 


1.  REPORT  DATE  (DD-MM-YYYY) 

2.  REPORT  TYPE 

September  2015 

Reprint 

4.  TITLE  AND  SUBTITLE 

On  Finsler  Geometry  and  Applications  in  Mechanics:  Review  and  New 
Perspectives 


6.  AUTHOR(S) 

JD  Clayton 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

US  Army  Research  Laboratory 

ATTN:  RDRL-WMP-C 

Aberdeen  Proving  Ground,  MD  21005-5069 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


11 


Hindawi  Publishing  Corporation 
Advances  in  Mathematical  Physics 
Volume  20 15,  Article  ID  828475, 11  pages 
http://dx.doi.org/10.1155/2015/828475 


<S> 

Hindawi 


Review  Article 

On  Finsler  Geometry  and  Applications  in  Mechanics: 
Review  and  New  Perspectives 


J,  D.  Clayton^’^ 

^Impact  Physics,  US  ARL,  Aberdeen,  MD  21005-5066,  USA 

^A.  James  Clark  School  of  Engineering  (Adjunct  Faculty),  University  of  Maryland,  College  Park,  MD  20742,  USA 
Correspondence  should  be  addressed  to  J.  D.  Clayton;  john.d.claytonl.civ@mail.mil 
Received  21  November  2014;  Accepted  18  January  2015 
Academic  Editor:  Mahouton  N.  Hounkonnou 

Copyright  ©  2015  J.  D.  Clayton.  This  is  an  open  access  article  distributed  under  the  Creative  Commons  Attribution  License,  which 
permits  unrestricted  use,  distribution,  and  reproduction  in  any  medium,  provided  the  original  work  is  properly  cited. 

In  Finsler  geometry,  each  point  of  a  base  manifold  can  be  endowed  with  coordinates  describing  its  position  as  well  as  a  set  of  one  or 
more  vectors  describing  directions,  for  example.  The  associated  metric  tensor  may  generally  depend  on  direction  as  well  as  position, 
and  a  number  of  connections  emerge  associated  with  various  covariant  derivatives  involving  affine  and  nonlinear  coefficients. 
Finsler  geometry  encompasses  Riemannian,  Euclidean,  and  Minkowskian  geometries  as  special  cases,  and  thus  it  affords  great 
generality  for  describing  a  number  of  phenomena  in  physics.  Here,  descriptions  of  finite  deformation  of  continuous  media  are  of 
primary  focus.  After  a  review  of  necessary  mathematical  definitions  and  derivations,  prior  work  involving  application  of  Finsler 
geometry  in  continuum  mechanics  of  solids  is  reviewed.  A  new  theoretical  description  of  continua  with  micro  structure  is  then 
outlined,  merging  concepts  from  Finsler  geometry  and  phase  field  theories  of  materials  science. 


1.  Introduction 

Mechanical  behavior  of  homogeneous  isotropic  elastic  solids 
can  be  described  by  constitutive  models  that  depend  only 
on  local  deformation,  for  example,  some  metric  or  strain 
tensor  that  may  generally  vary  with  position  in  a  body. 
Materials  with  microstructure  require  more  elaborate  consti¬ 
tutive  models,  for  example,  describing  lattice  orientation  in 
anisotropic  crystals,  dislocation  mechanisms  in  elastic-plastic 
crystals,  or  cracks  or  voids  in  damaged  brittle  or  ductile 
solids.  In  conventional  continuum  mechanics  approaches, 
such  models  typically  assign  one  or  more  time-  and  position- 
dependent  vector(s)  or  higher-order  tensor(s),  in  addition 
to  total  deformation  or  strain,  that  describe  physical  mech¬ 
anisms  associated  with  evolving  internal  structure. 

Mathematically,  in  classical  continuum  physics  [1-3], 
geometric  field  variables  describing  behavior  of  a  simply  con¬ 
nected  region  of  a  body  depend  fundamentally  only  on  ref¬ 
erential  and  spatial  coordinate  charts  {X^}  and  {x“}  {A,  a  = 
1, 2, . . .  n)  relatedby  adiffeomorphismx  =  (jp(X,t),  with  x  and 


X  denoting  corresponding  points  on  the  spatial  and  material 
manifolds  covered  by  corresponding  chart(s)  and  t  denoting 
time.  State  variables  entering  response  functions  depend 
ultimately  only  on  material  points  and  relative  changes  in 
their  position  (e.g.,  deformation  gradients  of  first  order 
and  possibly  higher  orders  for  strain  gradient-type  models 
[4]).  Geometric  objects  such  as  metric  tensors,  connection 
coefficients,  curvature  tensors,  and  anholonomic  objects  [5] 
also  depend  ultimately  only  on  position.  This  is  true  in 
conventional  nonlinear  elasticity  and  plasticity  theories  [1,6], 
as  well  as  geometric  theories  incorporating  torsion  and/or 
curvature  tensors  associated  with  crystal  defects,  for  example 
[7-15].  In  these  classical  theories,  the  metric  tensor  is  always 
Riemannian  (i.e.,  essentially  dependent  only  upon  x  or  X 
in  the  spatial  or  material  setting),  meaning  the  length  of  a 
differential  line  element  depends  only  on  position;  however, 
torsion,  curvature,  and/or  covariant  derivatives  of  the  metric 
need  not  always  vanish  if  the  material  contains  various 
kinds  of  defects  (non-Euclidean  geometry).  Connections  are 
linear  (i.e.,  affine).  Gauge  field  descriptions  in  the  context  of 
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Riemannian  metrics  and  affine  connections  include  [16, 17]. 
Relevant  references  in  geometry  and  mathematical  physics 
include  [18-26],  in  addition  to  those  already  mentioned. 
Finite  deformation  director  theories  of  micropolar  type  are 
addressed  in  the  context  of  Riemannian  (as  opposed  to 
Finslerian)  metrics  in  [1,  27]. 

Finsler  geometry,  first  attributed  to  Finsler  in  1918  [28], 
is  more  general  than  Riemannian  geometry  in  the  sense  that 
the  fundamental  (metric)  tensor  generally  may  depend  on 
additional  independent  variables,  labeled  here  as  y  and  Y 
in  spatial  and  material  configurations,  with  corresponding 
generalized  coordinates  [y“}  and  [Y^j.  Formal  definitions 
will  be  given  later  in  this  paper;  for  the  present  immediate 
discussion,  it  suffices  to  mention  that  each  point  can  be  con¬ 
sidered  endowed  with  additional  degrees-of-freedom  beyond 
X  or  X  and  that  transformation  laws  among  coordinates  as 
well  as  connection  coefficients  (i.e.,  covariant  differentials) 
generally  depend  on  y  or  y  as  well  as  x  or  X.  Relevant 
references  in  mathematics  include  [29-32] .  For  descriptions 
of  mechanics  of  solids,  additional  degrees-of-freedom  can  be 
associated  with  evolving  features  of  the  microstructure  of  the 
material,  though  more  general  physical  interpretations  are 
possible. 

The  use  of  Finsler  geometry  to  describe  continuum 
mechanical  behavior  of  solids  was  perhaps  first  noted  by 
Kroner  in  1968  [33]  and  Eringen  in  1971  [3],  the  latter 
reference  incorporating  some  basic  identities  and  definitions 
derived  primarily  by  Cartan  [34],  though  neither  devel¬ 
oped  a  Finsler-based  framework  more  specifically  directed 
towards  mechanics  of  continua.  The  first  theory  of  Finsler 
geometry  applied  to  continuum  mechanics  of  solids  with 
microstructure  appears  to  be  the  purely  kinematic  theory  of 
Ikeda  [35],  in  a  generalization  of  Cosserat-type  kinematics 
whereby  additional  degrees-of-freedom  are  director  vectors 
linked  to  structure.  This  theory  was  essentially  extended  by 
Bejancu  [30]  to  distinguish  among  horizontal  and  vertical 
distributions  of  the  fiber  bundle  of  a  deforming  pseudo- 
Finslerian  total  space.  More  complete  theories  incorporating 
a  Lagrangian  functional  (leading  to  physical  balance  or 
conservation  laws)  and  couched  in  terms  of  Finsler  geom¬ 
etry  were  developed  by  Saczuk,  Stumpf,  and  colleagues  for 
describing  solids  undergoing  inelastic  deformation  mecha¬ 
nisms  associated  with  plasticity  and/or  damage  [36-40].  To 
the  authors  knowledge,  solution  of  a  boundary  value  problem 
in  solid  mechanics  using  Finsler  geometric  theory  has  only 
been  reported  once,  in  [38].  Finsler  geometry  has  been 
analogously  used  to  generalize  fundamental  descriptions 
in  other  disciplines  of  physics  such  as  electromagnetism, 
quantum  theory,  and  gravitation  [30,  41-43]. 

This  paper  is  organized  as  follows.  In  Section  2,  requisite 
mathematical  background  on  Finsler  geometry  (sometimes 
called  Riemann-Finsler  geometry  [31])  is  summarized.  In 
Section  3,  the  aforementioned  theories  from  continuum 
physics  of  solids  [30,  35-38,  40]  are  reviewed  and  compared. 
In  Section  4,  aspects  of  a  new  theory,  with  a  primary  intention 
of  description  of  structural  transformation  processes  in  real 
materials,  are  proposed  and  evaluated.  Conclusions  follow  in 
Section  5. 


2,  Finsler  Geometry;  Background 

Notation  used  in  the  present  section  applies  to  a  referential 
description,  that  is,  the  initial  state;  analogous  formulae  apply 
for  a  spatial  description,  that  is,  a  deformed  body. 


2.1.  Coordinates  and  Fundamentals.  Denote  by  M  an  n- 
dimensional  C°°  manifold.  Each  element  (of  support)  of  M  is 
of  the  form  (X,  Y),  where  X  €  M  and  Y  e  TM,  with  TM  the 
tangent  bundle  of  M.  A  Finsler  structure  of  M  is  a  function 
L  :  TM  — >  [0,  oo)  with  the  following  three  properties  [31]: 

(i)  The  fundamental  function  L  is  C°°  on  TM  \  0; 

(ii)  L{X,  XY)  -  XL{X,  T)  VA  >  0  (i.e.,  L  is  homogeneous 
of  degree  one  in  T); 

(iii)  the  fundamental  tensor  =  {\l2)d^{l})ldY^dY^ 
is  positive  definite  at  every  point  of  TM  \  0. 

Restriction  of  L  to  a  particular  tangent  space  TxM  gives  rise 
to  a  (local)  Minkowski  norm 

l\Y)  =  Gj^{Y)Y^Y^,  (1) 


which  follows  from  Euler’s  theorem  and  the  identity 


Ld^ 

gyAgyiJ 


(2) 


Specifically  letting  Y^  — >  dX^,  the  length  of  a  differential 
line  element  at  X  depends  in  general  on  both  X  and  Y  as 

|dX  (X,  y)|  =  VdX  •  dX  =  [G^  (X,  Y)  dX^dX^]  .  (3) 

A  Finsler  manifold  (M,  F)  reduces  to  a  Minkowskian  man¬ 
ifold  when  L  does  not  depend  on  X  and  to  a  Riemannian 
manifold  when  L  does  not  depend  on  Y.  In  the  latter 
case,  a  Riemannian  metric  tensor  is  G^g  (X)dX^  (gi  dX^. 
Cartan’s  tensor,  with  the  following  fully  symmetric  covariant 
components,  is  defined  for  use  later: 

IdG^  1  (4) 

2  ayc  AdY^dY^dY^' 

Consider  now  a  coordinate  transformation  to  another 
chart  on  M;  for  example. 


X^  =  x^(x^x^...,x"), 


(5) 


From  the  chain  rule,  holonomic  basis  vectors  on  TM  then 
transform  as  [30, 31] 


a 

ax^ 


+  _a^^~c_a_ 
ax^  ax^  ax^^axc  ay^’ 

a  _  ax^  a 

ay^  ”  ax^ay^' 


(6) 

(7) 
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2.2.  Connections  and  Differentiation.  ChristofFel  symbols  of 
the  second  kind  derived  from  the  symmetric  fundamental 
tensor  are 


Ybc  -  2*^  dX^  dX°  )  ■ 


(8) 


Lowering  and  raising  of  indices  are  enabled  via  Gab  and  its 
inverse  Nonlinear  connection  coefficients  on  TM  \  0  are 
defined  as 


(9) 


where  =  yg^Y^Y‘".  The  following  nonholonomic  bases 
are  then  introduced: 


been  shown  [45]  how  components  of  Cartans  connection  on 
a  Finsler  manifold  can  be  obtained  as  the  induced  connection 
of  an  enveloping  space  (with  torsion)  of  dimension  2n.  When 
a  Finsler  manifold  degenerates  to  a  locally  Minkowski  space 
(L  independent  of  X),  then  F^^  =  =  0.  Gradients  of 

bases  with  respect  to  the  Chern-Rund  connection  and  Cartan 
tensor  are 


S  _c  5 
S/SX'^gXB  ~  ABgxC’ 


^Sisx^^  -  (16) 


As  an  example  of  covariant  differentiation  on  a  Finsler 
manifold  with  Chern-Rund  connection  V,  consider  a  (J) 
tensor  field  T  =  Tg{dldX^)  iS>  dX^  on  the  manifold  TM  \  0. 
The  covariant  differential  of  T(X,  Y)  is 


5 

dx^ 


d 


dx^.  (10) 


It  can  be  shown  that  unlike  (6),  these  nonholonomic 
bases  obey  simple  transformation  laws  like  (7).  The  set 
(d/dX^,  djdY^}  serves  as  a  convenient  local  basis  for  T{TM  \ 
0);  its  dual  set  [dX^,  dT^}  applies  for  the  cotangent  bundle 
T*{TM  \  0).  A  natural  Riemannian  metric  can  then  be 
introduced,  called  a  Sasaki  metric  [31]: 


G  (X,  y)  =  G^gdx^  g>  dx^  +  G^dy^  g>  dy^.  (ii) 


The  horizontal  subspace  spanned  by  [d/dX^j  is  orthogonal 
to  the  vertical  subspace  spanned  by  {djdY^}  with  respect  to 
this  metric.  Covariant  derivative  V,  or  collectively  connection 
1- forms  co^,  define  a  linear  connection  on  pulled-back  bundle 
n*TM  over  TM  \  0.  Letting  v  denote  an  arbitrary  direction, 

V.  <  (^)  ^^dX'^  =  -"b  (v)  dx\  (12) 


A  number  of  linear  connections  have  been  introduced  in 
the  Finsler  literature  [30,  31].  The  Chern-Rund  connection 
[29,  44]  is  used  most  frequently  in  applications  related  to 
the  present  paper.  It  is  a  unique  linear  connection  on  n*TM 
characterized  by  the  structural  equations  [31] 


d  (dX"^)  -  dX^  A  a)g  =  0, 


^i^AB  ^BC^A  ^AC^B  ~  ^^ABC^^  ■ 


(13) 


The  first  structure  equation  implies  torsion  freeness  and 
results  in 


A  T'-'^  J 


_  ■p'-A 

^BC  ~  ^CB' 


(14) 


The  second  leads  to  the  connection  coefficients 


^BC 


SG 


BD 


SG, 


CD 


SG 


BC 


dxc  dx^  dx^  J 


). 


(15) 


When  a  Finsler  manifold  degenerates  to  a  Riemannian 
manifold,  =  0  and  F^(^  =  y^^.  Cartans  connection  1- forms 
are  defined  by  to^  +  C^^SY^  where  o)^  correspond  to  (14); 
its  coordinate  formulae  and  properties  are  listed  in  [3].  It  has 


{VT)^  =  dT^  +  T^4- 


ri-iA  C 
^C^B 


-{VsisxcT)^dX^  +  {Vs/3YcT)^SY^ 


(STl 

Vdxc 


ri-iDr'A  rT-iAr'D 

^B  ^CD  ~  ^D^CB 


dX^  -t 


(djl 

I  dy*^ 


dy*^. 

(17) 


Notations  (Oia  ‘‘nd  (Oja  denote  respective  horizontal  and 
vertical  covariant  derivatives  with  respect  to  V. 


2.3.  Geometric  Quantities  and  Identities.  Focusing  again  on 
the  Chern-Rund  connection  V,  curvature  2-forms  are 


,-^A  AN  C  . 

Llg  -  a  Ico^  1  -  o)^  A  o). 


^codX^  A  dX°  +  A  SY^  (ig) 


+  ^Qbcd^y''^SY^ 


with  d(-)  the  exterior  derivative  and  A  the  wedge  product  (no 
factor  of  1/2).  HH-,  HV-,  and  VV-curvature  tensors  of  the 
Chern-Rund  connection  have  respective  components 


R 


ST 


BD 


ST, 


BCD 


‘BCD 


dxc  dx^ 

_^£bC  ( 

ay^’ 


BC  ,  •pA  p£  _  pA  p£ 

^EC^BD  ^ED^BC^ 


<BCD 


=  0. 


(19) 


VV-curvature  vanishes,  HV-curvature  obeys  P^^d  ~  ^ci 
and  a  Bianchi  identity  for  HH-curvature  is 


I  pA  ,  pA  _  ^ 
^BCD  ^CDB  ^DBC  ~ 


(20) 


When  a  Finsler  manifold  degenerates  to  a  Riemannian 
manifold,  then  R^^d  become  the  components  of  the  usual 
curvature  tensor  of  Riemannian  geometry  constructed  from 
and  P^cd  ~  b-  AH  curvatures  vanish  in  locally 
Minkowski  spaces.  It  is  not  always  possible  to  embed  a  Finsler 
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space  in  a  Riemannian  space  without  torsion,  but  it  is  possible 
to  determine  the  metric  and  torsion  tensors  of  a  space  of 
dimension  2n  -  1  in  such  a  way  that  any  «- dimensional 
Finsler  space  is  a  nonholonomic  subspace  of  such  a  space  with 
torsion  [46]. 

Nonholonomicity  (i.e.,  nonintegrability)  of  the  horizontal 
distribution  is  measured  by  [32] 


S  S 


SX'^’SXB 


SNi  A_  ™ 


sx^  sx^  )  ayc 


where  [■,  ■]  is  the  Lie  bracket  and  can  be  interpreted  as 
components  of  a  torsion  tensor  [30].  For  the  Chern-Rund 
connection  [31], 


^'-AB  ~  ^DAB^  ■ 


(22) 


Since  Lie  bracket  (21)  is  strictly  vertical,  the  horizontal 
distribution  spanned  by  [5/5X^}  is  not  involutive  [31]. 


3.  Applications  in  Solid  Mechanics:  1973-2003 


with  =  B^(X).  Differentiating  the  second  of  (25), 

+  ^bK):a)  (26) 

where  (O-a  denotes  the  total  covariant  derivative  [20,  26]. 
Differential  line  and  director  elements  can  be  related  by 


Let  Gj^{X)  denote  the  metric  in  the  reference  con¬ 
figuration,  such  that  dS^  =  dX^GA^dX^  is  a  measure 
of  length  (Gab  “  ^ab  for  Cartesian  coordinates  [X^}). 
Fundamental  tensors  in  the  spatial  frame  describing  strains 
of  the  continuum  and  directors  are 

'^ab  -  ^ABh  ’  '^ab  “  ^{a)a^AB^{a)b- 

Let  F^p  and  F^^  denote  coefficients  of  linear  connections 
associated  with  continuum  and  director  fields,  related  by 


3.1.  Early  Director  Theory.  The  first  application  of  Finsler 
geometry  to  finite  deformation  continuum  mechanics  is 
credited  to  Ikeda  [35] ,  who  developed  a  director  theory  in  the 
context  of  (pseudo-)  Finslerian  manifolds.  A  slightly  earlier 
work  [47]  considered  a  generalized  space  (not  necessarily 
Finslerian)  comprised  of  finitely  deforming  physical  and 
geometrical  fields.  Paper  [35]  is  focused  on  kinematics  and 
geometry:  descriptions  of  deformations  of  the  continuum 
and  the  director  vector  fields  and  their  possible  interactions, 
metric  tensors  (i.e.,  fundamental  tensors),  and  gradients  of 
motions.  Covariant  differentials  are  defined  that  can  be  used 
in  field  theories  of  Finsler  space  [41,  48].  Essential  concepts 
from  [35]  are  reviewed  and  analyzed  next. 

Let  M  denote  a  pseudo-Finsler  manifold  in  the  sense  of 
[35],  representative  of  a  material  body  with  microstructure. 
Let  X  6  M  denote  a  material  point  with  coordinate  chart 
[X^l,  (A  =  1, 2, 3)  covering  the  body  in  its  undeformed  state. 
A  set  of  director  vectors  is  attached  to  each  X,  where 
(a  =  1, 2, . . . ,  p).  In  component  form,  directors  are  written 
[D^jj.  In  the  context  of  notation  in  Section  2,  M  is  similar 
to  a  Finsler  manifold  with  y^  ^  D^j,  though  the  director 
theory  involves  more  degrees-of-freedom  when  p  >  l,andno 
fundamental  function  is  necessarily  introduced.  Let  x  denote 
the  spatial  location  in  a  deformed  body  of  a  point  initially 
at  X,  and  let  d(„)  denote  a  deformed  director  vector.  The 
deformation  process  is  described  by  (a  =  1,2, 3) 

(X^)  ,  dl,  =  dl,  [Df^,  (X^)]  =  dl,  (X")  . 

(23) 


(24) 


(25) 


The  deformation  gradient  and  its  inverse  are 

pa  -lA  ^  ^ 

^  ax^’  “  ■ 

The  following  decoupled  transformations  are  posited: 


Tbc  - 


(29) 


where  D^jD^“'  =  dg  and  D^jD^'  =  5^.  The  covariant 
differential  of  a  referential  vector  field  V,  where  locally 
V(X)  e  T^M,  with  respect  to  this  connection  is 


[vv)'^  =  dv^  +  fX'^dx^;  (30) 


the  covariant  differential  of  a  spatial  vector  field  v,  where 
locally  v{x)  e  T^M,  is  defined  as 

(Vy)  =  dv'  4-  F^^y'^dx*'  +  F|,*''’y‘^dx*',  (31) 


with,  for  example,  the  differential  of  y“  (x,  d(„) )  given  by  dy“  = 
(ay“/ax^)dx^  +  (ay“/a(i(„))d<i(^).  Euler-Schouten  tensors  are 


j_r{a)a  _  jnU  -p-iA 

^bc  ~  ^A^(a)c:b> 


^{a)a 

^bc 


“  ^{a)A^c:b  * 


(32) 


Ikeda  [35]  implies  that  fundamental  variables  entering  a 
field  theory  for  directed  media  should  include  the  set 
(F,  B,E,  H,  K).  Given  the  fields  in  (23),  the  kinematic- 

— Ct 

geometric  theory  is  fuUy  determined  once  F^^  are  defined. 
The  latter  coefficients  can  be  related  to  defect  content  in 
a  crystal.  For  example,  setting  F^^  =  0  results  in  distant 
parallelism  associated  with  dislocation  theory  [7],  in  which 
case  (29)  gives  the  negative  of  the  wryness  tensor;  more 
general  theory  is  needed,  however,  to  represent  disclination 
defects  [49]  or  other  general  sources  of  incompatibility  [6, 
50]  requiring  a  nonvanishing  curvature  tensor.  The  directors 
themselves  can  be  related  to  lattice  directions,  slip  vectors  in 
crystal  plasticity  [51,  52],  preferred  directions  for  twinning 
[6,  53],  or  planes  intrinsically  prone  to  cleavage  fracture  [54], 
for  example.  Applications  of  multifield  theory  [47]  towards 
multiscale  descriptions  of  crystal  plasticity  [50,  55]  are  also 
possible. 


Advances  in  Mathematical  Physics 


5 


3.2.  Kinematics  and  Gauge  Theory  on  a  Fiber  Bundle.  The 
second  known  application  of  Finsler  geometry  towards  finite 
deformation  of  solid  bodies  appears  in  Chapter  8  of  the  book 
of  Bejancu  [30].  Content  in  [30]  extends  and  formalizes  the 
description  of  Ikeda  [35]  using  concepts  of  tensor  calculus  on 
the  fiber  bundle  of  a  (generalized  pseudo-)  Finsler  manifold. 
Geometric  quantities  appropriate  for  use  in  gauge-invariant 
Lagrangian  functions  are  derived.  Relevant  features  of  the 
theory  in  [30]  are  reviewed  and  analyzed  in  what  follows  next. 

Define  C  =  (-Z,  n,  M,  U)  as  a  fiber  bundle  of  total  space 
Z,  where  n  :  Z  M  is  the  projection  to  base  manifold 
M  and  U  is  the  fiber.  Dimensions  of  M  and  U  are  n  {n  -  3 
for  a  solid  volume)  and  p,  respectively;  the  dimension  of  Z  is 
r  =  n  +  p.  Coordinates  on  Z  are  [X^,D“},  where  X  €  M 
is  a  point  on  the  base  body  in  its  reference  configuration, 
and  D  is  a  director  of  dimension  p  that  essentially  replaces 
multiple  directors  of  dimension  3  considered  in  Section  3.1. 
The  natural  basis  on  Z  is  the  field  of  frames  {d/dX^,  djdD'^}. 
Let  N'^iX,  D)  denote  nonlinear  connection  coefficients  on  Z, 
and  introduce  the  nonholonomic  bases 


5 

SX^ 


d 

90“’ 


5D“  =  dD‘"  +  N'ldX^. 


(33) 


Unlike  {djdX^},  these  nonholonomic  bases  obey  simple 
transformation  laws;  [5/5X^,  9/9D“}  serves  as  a  convenient 
local  basis  for  TZ  adapted  to  a  decomposition  into  horizontal 
and  vertical  distributions: 


TZ^  HZ®  VZ. 


(34) 


Generalization  to  higher-order  tensor  fields  is  given  in  [30]. 
In  particular,  the  following  coefficients  are  assigned  to  the  so- 
called  gauge  H-connection  on  Z: 

pA  _  }_(^AD  (  ^ 

2  \  8X^  SX^  J  ’ 


pa  dK 

dD^  ’ 

X^y=zG  I  ^  + 


C  -  0, 


9G^3  dG 


(39) 


Ty 


doy  dD^  dD^ 


Gomparing  with  the  formal  theory  of  Finsler  geometry 
outlined  in  Section  2,  coefficients  F^(^  are  analogous  to  those 
of  the  Ghern-Rund  connection,  and  (36)  is  analogous  to 
(16).  The  generalized  pseudo-Finslerian  description  of  [30] 
reduces  to  Finsler  geometry  of  [31]  when  p  -  n  and  a 
fundamental  function  L  exists  from  which  metric  tensors  and 
nonlinear  connection  coefficients  can  be  derived. 

Let  Q^{X,  D)  denote  a  set  of  differentiable  state  variables, 
where  K  =  1,2, ...,r.  A  Lagrangian  function  2  of  the 
following  form  is  considered  on  Z: 


2(X,D)  =  2 


“AB 


^  '  8X^'  9D“ 


(40) 


Let  D  be  a  compact  domain  of  Z,  and  define  the  functional 
(action  integral) 


7  (D)  =  [  2  (X,  D)  dx'  •  •  •  A  dX"  A  dD*  •  •  •  A  dD^, 

Jn 


(41) 


2=  [|det(G^5)  II  det(G„^)|]'''2. 


Dual  set  [dX^,  5D“}  is  a  local  basis  on  r*Z.  A  fundamental 
tensor  (i.e.,  metric)  for  the  undeformed  state  is 


Euler-Lagrange  equations  referred  to  the  reference  configu¬ 
ration  follow  from  the  variational  principle  8I  =  0: 


92 

9 

92 

9 

92 

G  (X,  D)  =  Gab  (X,  D)  dX^  (g>  dX^ 

9Q^ 

dx^ 

9(9Q^/9XA) 

9D“ 

9  {dQ^ldD‘^) 

■G„8  (X,D)5D“g>dD^. 


(35) 


Let  V  denote  co variant  differentiation  with  respect  to  a 
connection  on  the  vector  bundles  HZ  and  VZ,  with 


5 


—  _ 

9 

-c^ 

^dxc’ 

^S/SX^  9D“ 

Aa  Qjjj}  ’ 

B  ^ 

„  9 

r  a 

Xaji^Dy' 

^a/3D“  ■ 


Gonsider  a  horizontal  vector  field  V  =  V^{dldX^)  and  a 
vertical  vector  field  W  =  W“(9/9D“).  Horizontal  and  vertical 
covariant  derivatives  are  defined  as 


A  AC  a  SW“  oc  B 

A  SV^  AC  a  „  s 


(38) 


(42) 

These  can  be  rewritten  as  invariant  conservation  laws  involv¬ 
ing  horizontal  and  vertical  covariant  derivatives  with  respect 
to  the  gauge-H  connection  [30]. 

Let  ('  -  {Z' ,n' ,m' ,u')  be  the  deformed  image  of 
fiber  bundle  C>  representative  of  deformed  geometry  of  the 
body,  for  example.  Dimensions  of  M'  and  U'  are  n  and  p, 
respectively;  the  dimension  of  Z'  isr  -  n  +  p.  Goordinates  on 
Z  are  [x“,  <i“  },  where  x  6  M  is  a  point  on  the  base  body  in 
its  current  configuration,  and  d  is  a  director  of  dimension  p. 

The  natural  basis  on  Z'  is  the  field  of  frames  [9/9x“,  djdd'^  }. 
Let  ]V“  (x,  d)  denote  nonlinear  connection  coefficients  on  Z', 
and  introduce  the  nonholonomic  bases 

^  ^  -  ]V“'  dd“'  =  dd“'  +  Kdx^  (43) 

5x“  9x“  “  9d“' 

These  nonholonomic  bases  obey  simple  transformation  laws; 

[d/dx“,9/9d“  }  serves  as  a  convenient  local  basis  for  TZ' 
where 

TZ'  =  HZ'  ©  VZ'. 


(44) 
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Dual  set  {dx‘‘,8cr  }  is  a  local  basis  on  ^*Z^  Deformation 
of  (Z,  M)  to  {Z' ,  M')  is  dictated  by  diffeomorphisms  in  local 
coordinates: 

(X^)  ,  d“'  =  (X)  D^.  (45) 


Let  the  usual  (horizontal)  deformation  gradient  and  its 
inverse  have  components 


Pa  = 


dx‘‘ 

dX^’ 


dX^ 

dx" 


(46) 


It  follows  from  (45)  that,  similar  to  (6)  and  (7), 


d 

dX^ 


db; 

dX^ 


13' 


d 

dd^'' 


(47) 


Nonlinear  connection  coefficients  on  Z  and  Z'  can  be  related 
by 


K  Pa  =  N-b;  - 


'aB“  ' 


D'’ 


(48) 


of  [30]  and  Section  3.2  to  describe  mechanics  of  solids  with 
microstructure  undergoing  finite  elastic-plastic  or  elastic- 
damage  deformations.  Key  new  contributions  of  these  works 
include  definitions  of  total  deformation  gradients  consist¬ 
ing  of  horizontal  and  vertical  components  and  Lagrangian 
functions  with  corresponding  energy  functionals  dependent 
on  total  deformations  and  possibly  other  state  variables. 
Constitutive  relations  and  balance  equations  are  then  derived 
from  variations  of  such  functionals.  Essential  details  are 
compared  and  analyzed  in  the  following  discussion.  Notation 
usually  follows  that  of  Section  3.2,  with  a  few  generalizations 
defined  as  they  appear. 

Define  C  =  (-Z,  n,  M,  U)  as  a  fiber  bundle  of  total  space 
Z,  where  n  :  Z  —>■  M  is  the  projection  to  base  manifold 
M  and  U  is  the  fiber.  Dimensions  of  M  and  U  are  3  and  p, 
respectively;  the  dimension  of  Z  is  r  =  3  +  p.  Coordinates 
on  Z  are  (X^,  D“},  where  X  6  M  is  a  point  on  the  base  body 
in  its  reference  configuration,  and  D  is  a  vector  of  dimension 
p.  Let  ('  =  {Z' ,n  ,M' ,U')  be  the  deformed  image  of  fiber 
bundle  C;  dimensions  of  M'  and  U'  are  3  and  p,  respectively; 
the  dimension  of  Z^  is  r  =  3  +  p.  Coordinates  on  Z'  are 

{x‘‘,d“  },  where  x  6  M  is  a  point  on  the  base  body  in 
its  current  configuration,  and  d  is  a  vector  of  dimension  p. 
Tangent  bundles  can  be  expressed  as  direct  sums  of  horizontal 
and  vertical  distributions: 


Metric  tensor  components  on  Z  and  Z'  can  be  related  by 

GahPlpB  =  G^,  Bf  =  G,,.  (49) 

Linear  connection  coefficients  on  Z  and  Z'  can  be  related  by 


TZ  =  HZ®  VZ,  TZ'  =  HZ'  ©  VZ' .  (51) 

Deformation  of  Z  to  Z'  is  locally  represented  by  the  smooth 
and  invertible  coordinate  transformations 


T-iU  jjC 

^bc^B^C  ~  ^BC^A 


(x"^,  D“)  ,  d“'  =  B^'d^,  (52) 


'  a'  u  ' 

a  r^b  r^cc  B  r~iCi 

<3a'hpABl3  =CfiApB^ 


(50) 


where  in  general  the  director  deformation  map  [38] 


B“  =B“  (X,D). 


(53) 


Total  deformation  gradient  S  :  PZ  TZ'  is  defined  as 


Using  the  above  transformations,  a  complete  gauge  Id- 
connection  can  be  obtained  for  Z'  from  reference  quantities 
on  Z  if  the  deformation  functions  in  (45)  are  known.  A 
Lagrangian  can  then  be  constructed  analogously  to  (40),  and 
Euler-Lagrange  equations  for  the  current  configuration  of 
the  body  can  be  derived  from  a  variational  principle  where 
the  action  integral  is  taken  over  the  deformed  space.  In 
application  of  pseudo-Finslerian  fiber  bundle  theory  similar 
to  that  outlined  above,  Fu  et  al.  [37]  associate  d  with  the 
director  of  an  oriented  area  element  that  may  be  degraded 
in  strength  due  to  damage  processes  such  as  fracture  or  void 
growth  in  the  material.  See  also  related  work  in  [39]. 


S  -  F  +  F  -  ®  ^  ® 


In  component  form,  its  horizontal  and  vertical  parts  are 


U*  a  ,  rB  cC  cfl  c 

Pa  ~  ^\A  ~  ^B^  ’ 

pl3'  _  -I-  S^S‘'S‘^x‘'\  s^' 


(55) 


3.3.  Recent  Theories  in  Damage  Mechanics  and  Finite  Plastic¬ 
ity.  Saczuk  et  al.  [36, 38, 40]  adapted  a  generalized  version  of 
pseudo-Finsler  geometry  similar  to  the  fiber  bundle  approach 


To  eliminate  further  excessive  use  of  Kronecker  deltas,  let 
p  -  3,  -  5^D“,  and  d‘‘  -  d'^id'" .  Introducing 

a  fundamental  Finsler  function  L{X,D)  with  properties 
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described  in  Section  2.1  {Y  —>■  D),  the  following  definitions 
hold  [38] : 


2  dD^dD^  ’ 

A  _  ^  dG^c  , 

Ybc  -2^  axc  dx^  dX°  )  ’ 


6 

SX^ 


A  1  3G^  d  ABC 

Nr  =  -z—s,  G^  =  YnrD  D^; 

^  2dD^ 

2  dDC  AdD^dD^dD^’ 
pA  _  ^(-.AD  /  fiG^p  dG(^p  ^  dG^(^  \ 

2  V  dxc  dx^  dx°  j’ 

=  ^  -nI^>  SD^  =  dD^  +  Ngdx’^. 
dX^  ^ 


(56) 


Notice  that  correspond  to  the  Chern-Rund  connection 
and  G^c  to  the  Cartan  tensor.  In  [36],  D  is  presumed 

t  r 

stationary  (B^  ^  d^  )  and  is  associated  with  a  residual 

plastic  disturbance,  G^^  is  remarked  to  be  associated  with 
dislocation  density,  and  the  HV-curvature  tensor  of  (e.g., 
Pg(jD  (^^))  remarked  to  be  associated  with  disclination 
density.  In  [38],  B  is  associated  with  lattice  distortion  and 
with  residual  strain  energy  density  of  the  dislocation  density 
[6,  56].  In  [40],  a  fundamental  function  and  connection 
coefficients  are  not  defined  explicitly,  leaving  the  theory  open 
to  generalization.  Note  also  that  (52)  is  more  general  than 
(45).  When  the  latter  holds  and  P^^  =  0,  then  g  ^  F>  the 
usual  deformation  gradient  of  continuum  mechanics. 

Restricting  attention  to  the  time-independent  case,  a 
Lagrangian  is  posited  of  the  form 


S(X,D) 

=  S  [X,  D,  X  (X,  D) ,  g  (X,  D) ,  Q  (X,  D) ,  VQ  (X,  D)] , 

(57) 


where  Q  is  a  generic  vector  of  state  variables  and  V  denotes  a 
generic  gradient  that  may  include  partial,  horizontal,  and/or 
vertical  covariant  derivatives  in  the  reference  configuration 
as  physically  and  mathematically  appropriate.  Let  fl  be  a 
compact  domain  of  Z,  and  define 


P  (fl)  =  [  a  (X,  D)  dX*  •  •  •  A  dX"  A  do'  •  •  •  A  dD^, 

Jn 


|det(G^)| 


-,1/2 


(58) 


a. 


Euler-Lagrange  equations  referred  to  the  reference  configura¬ 
tion  follow  from  the  variational  principle  dl  -  0  analogously 


to  (42).  Conjugate  forces  to  variations  in  kinematic  and  state 
variables  can  be  defined  as  derivatives  of  2  with  respect 
to  these  variables.  Time  dependence,  dissipation,  first  and 
second  laws  of  thermodynamics,  and  temperature  effects  are 
also  considered  in  [38, 40] ;  details  are  beyond  the  scope  of  this 
review.  In  the  only  known  application  of  Finsler  geometry  to 
solve  a  boundary  value  problem  in  the  context  of  mechanics 
of  solids  with  microstructure,  Stumpf  and  Saczuk  [38]  use 
the  theory  outlined  above  to  study  localization  of  plastic  slip 
in  a  bar  loaded  in  tension,  with  specifying  a  preferred 
material  direction  for  slip.  In  [40],  various  choices  of  Q  and 
its  gradient  are  considered  in  particular  energy  functions,  for 
example,  state  variables  associated  with  gradients  of  damage 
parameters  or  void  volume  fractions. 

4.  Towards  a  New  Theory  of  Structured  Media 

4.1.  Background  and  Scope.  In  Section  4  it  is  shown  how 
Finsler  geometry  can  be  applied  to  describe  physical  prob¬ 
lems  in  deformable  continua  with  evolving  microstructures 
in  a  manner  somewhat  analogous  to  the  phase  field  method. 
Phase  field  theory  [57]  encompasses  various  diffuse  interface 
models  wherein  the  boundary  between  two  (or  more)  phases 
or  states  of  material  is  distinguished  by  the  gradient  of  a  scalar 
field  called  an  order  parameter.  The  order  parameter,  denoted 
herein  by  rj,  typically  varies  continuously  between  values  of 
zero  and  unity  in  phases  one  and  two,  with  intermediate 
values  in  phase  boundaries.  Mathematically, 

tj  (X)  =  0  VX  6  phase  1, 

(X)  =  1  VX  e  phase  2,  (59) 

tj  (X)  6  (0, 1)  VX  6  interface. 

Physically,  phases  might  correspond  to  liquid  and  solid 
in  melting- solidification  problems,  austenite  and  martensite 
in  structure  transformations,  vacuum  and  intact  solid  in 
fracture  mechanics,  or  twin  and  parent  crystal  in  twinning 
descriptions.  Similarities  between  phase  field  theory  and 
gradient-type  theories  of  continuum  mechanics  are  described 
in  [58];  both  classes  of  theory  benefit  from  regularization 
associated  with  a  length  scale  dependence  of  solutions  that 
can  render  numerical  solutions  to  governing  equations  mesh 
independent.  Recent  phase  field  theories  incorporating  finite 
deformation  kinematics  include  [59]  for  martensitic  trans¬ 
formations,  [54]  for  fracture,  [60]  for  amorphization,  and 
[53]  for  twinning.  The  latter  (i.e.,  deformation  twinning)  is 
the  focus  of  a  more  specific  description  that  follows  later  in 
Section  4.3  of  this  paper,  though  it  is  anticipated  that  the 
Finsler-type  description  could  be  adapted  straightforwardly 
to  describe  other  deformation  physics. 

Deformation  twinning  involves  shearing  and  lattice  rota¬ 
tion/reflection  induced  by  mechanical  stress  in  a  solid  crystal. 
The  usual  elastic  driving  force  is  a  resolved  shear  stress  on 
the  habit  plane,  in  the  direction  of  twinning  shear.  Twinning 
can  be  reversible  or  irreversible  depending  on  material  and 
loading  protocol;  the  physics  of  deformation  twinning  is 
described  more  fully  in  [61]  and  Chapter  8  of  [6] .  Theory  that 
follows  in  Section  4.2  is  thought  to  be  the  first  to  recognize 


8 


Advances  in  Mathematical  Physics 


analogies  between  phase  field  theory  and  Finsler  geometry 
and  that  in  Section  4.3  the  first  to  apply  Finsler  geometric 
concepts  to  model  deformation  twinning.  Notation  follows 
that  of  Section  3.2  and  Section  3.3  with  possible  exceptions 
highlighted  as  they  appear. 

4.2.  Finsler  Geometry  and  Kinematics.  As  in  Section  3.3, 
define  (  -  (Z,  n,  M,  U)  as  a  fiber  bundle  of  total  space  Z 
with  n  :  Z  M.  Considered  is  a  three-dimensional  solid 
body  with  state  vector  D:  dimensions  of  M  and  U  are  3; 
the  dimension  of  Z  is  6.  Coordinates  on  Z  are 
{A  -  1, 2, 3),  where  X  6  M  is  a  point  on  the  base  body 
in  its  reference  configuration.  Let  ('  =  be 

the  deformed  image  of  fiber  bundle  C;  dimensions  of  M' 
and  U'  are  3,  and  that  of  Z'  is  6.  Coordinates  on  Z'  are 
where  x  6  M'  is  a  point  on  the  base  body  in 
its  current  configuration,  and  d  is  the  updated  state  vector. 
Tangent  bundles  are  expressed  as  direct  sums  of  horizontal 
and  vertical  distributions  as  in  (51): 

TZ  =  HZ  e  VZ,  TZ’  =  HZ’  ©  VZ’ .  (60) 


From  (63)  and  (65),  local  integrability  (null  curl)  conditions 
dD^jdX^  =  dD^jdX^  and  ddVdx*’  =  hold. 

As  in  Section  3,  Ng{X,  D)  denote  nonlinear  connection 
coefficients  on  Z,  and  the  Finsler-type  nonholonomic  bases 


d 


3  3 


8X^  dX^ 


SD'^  =  dD'^  +  NgdX^.  (67) 


Let  =  NgS‘^Sy  denote  nonlinear  connection  coefficients 
on  Z'  [38],  and  define  the  nonholonomic  spatial  bases  as 


dx“  dx“  “  dd^  ^ 


(68) 


Total  deformation  gradient  S  :  TZ  TZ'  is  defined  as 


In  component  form,  its  horizontal  and  vertical  parts  are,  as  in 
(55)  and  the  theory  of  [38], 


Deformation  of  Z  to  Z'  is  locally  represented  by  the  smooth 
and  invertible  coordinate  transformations 

x“  =  x“  (X^,  D^)  ,  d“  =  (61) 

where  the  director  deformation  function  is,  in  general, 

B“-Bl(X,D).  (62) 


a  dx  D  g  j 

Ta  ~  ^|A  “  ^  > 


pa  _  a  _  c 

-  ^\\A  -  Qj^A  ^  '-"AC°bOc  ^  • 


(70) 


As  in  Section 3.3,  a  fundamental  Finsler  function  L{X,D), 
homogeneous  of  degree  one  in  D,  is  introduced.  Then 


Introduce  the  smooth  scalar  order  parameter  field  rj  €  M 
as  in  (59),  where  rj  -  r]{X).  Here  D  is  identified  with  the 
reference  gradient  of  rj: 

D^{X)  =  d^^^.  (63) 


For  simplicity,  here  [X^j  is  taken  as  a  Cartesian  coordinate 
chart  on  M  such  that  X^  =  S^^Xg,  and  so  forth.  Similarly, 
[x“}  is  chosen  Cartesian.  Generalization  to  curvilinear  coor¬ 
dinates  is  straightforward  but  involves  additional  notation. 
Define  the  partial  deformation  gradient 


fUx,d)  = 


ax“  (x,D) 
ax^  ■ 


By  definition  and  then  from  the  chain  rule, 

^abdri  (X) 


[x(X,D),D]  =  d“ 


ax^ 


'  ax^  [dx^J 


(64) 


(65) 


^  1 

'-'AB  “  -  ■ 


2  dD^dD^  ’ 


_A  i  ^AD  (  ^Ggjj 


3G, 


CD 


dG 


BC 


dx^  dx^ 


^  2dD^ 


C 


1  ac/ 


1 


a^(L^) 


ABC 


2  aoc  AdD^dD^dD^’ 


pA  _  1  ^AD 

Tbc  -  2*^ 


dG 


BD 


dG, 


CD 


dG 


BC 


dxc  dx^  dx^ 


(71) 

(72) 

(73) 

(74) 

(75) 


Specifically  for  application  of  the  theory  in  the  context  of 
initially  homogeneous  single  crystals,  let 


AB 


(  dtj  \/  dtj  \ 

vax^y  vaxvj 


-,1/2 


-  Ga„-  dArdnnG’"^  -  Constant. 


•^AC'^BD' 


(76) 

(77) 


—  \A  ocjfc  c  t~\B 

b  3  d^D  , 

leading  to 

B“  (X,  D)  =  Fl^^  (X,  D)  8^'’8j,g.  (66) 


The  metric  tensor  with  Cartesian  components  is  iden¬ 
tified  with  the  gradient  energy  contribution  to  the  surface 
energy  term  in  phase  field  theory  [53] ,  as  will  be  made  explicit 
later  in  Section  4.3.  From  (76)  and  (77),  reference  geometry  C 
is  now  Minkowskian  since  the  fundamental  Finsler  function 
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L  =  L{D)  is  independent  of  X.  In  this  case,  (72)-(75)  and 
(67)-(68)  reduce  to 


^BC  ~  VbC  ~  '^BC  ~ 

S  _  d 
SX^  ~  dx^’ 

A  -  A 

’ 


Ci 

n.  =  0; 


=  dD^-, 


(78) 


=  dd‘‘. 


Horizontal  and  vertical  covariant  derivatives  in  (69)  reduce 
to  partial  derivatives  with  respect  to  X^  and  D^,  respectively, 
leading  to 


A 

dX^  dx‘‘ 
d 


dD^  9<i“ 
d 


)dD^ 


(79) 


=  —  «>  dX^  +  Vl—®  dD 


When  X  =  x(X),  then  vertical  deformation  components 
Pa  ~  ~  ^  more  general  version  of  (76)  applicable 

to  heterogeneous  material  properties,  Kab  =  Kab(X,  D)  and  is 
homogeneous  of  degree  zero  with  respect  to  D,  and  the  above 
simplifications  (i.e.,  vanishing  connection  coefficients)  need 
not  apply. 


4.3.  Governing  Equations:  Twinning  Application.  Consider  a 
crystal  with  a  single  potentially  active  twin  system.  Applying 
(59),  let  rj{X)  =  1  VX  6  twinned  domains,  rj{X)  =  0  VX  6 
parent  (original)  crystal  domains,  and  rj(X)  6  (0, 1)  VX  6  in 
twin  boundary  domains.  As  defined  in  [53],  let  denote  the 
twinning  eigenshear  (a  scalar  constant),  m  =  m-AdX^  denote 
the  unit  normal  to  the  habit  plane  (i.e.,  the  normal  covector 
to  the  twin  boundary),  and  s  =  5^(d/dX^)  the  direction  of 
twinning  shear.  In  the  context  of  the  geometric  framework 
of  Section  4.2,  and  with  simplifying  assumptions  (76)-(79) 
applied  throughout  the  present  application,  s  6  TM  and 
m  6  r*M  are  constant  fields  that  obey  the  orthonormality 
conditions 

(s, m)  =  ^  dX^^  ^  s^niA  =  0.  (80) 


Twinning  deformation  is  defined  as  the  (J)  tensor  field 


“  "  "  1  +  yo<ps«>m. 

P-B  [’I  (^)]  =Sb+  Vof  [’I  (^)] 


(81) 


Note  that  detH  =  1  +  yo<p(s>iii)  -  1.  Scalar  interpolation 
function  (p{rj)  6  [0, 1]  monotonically  increases  between  its 
endpoints  with  increasing  rj,  and  it  satisfies  (p(0)  -  0,  (jp(l)  = 
1,  and  (p'{0)  =  (p'{l)  =  0.  A  typical  example  is  the  cubic 
polynomial  [53,  59] 


The  elastic  lattice  deformation  is  the  two -point  tensor  A: 

A  =  ®  ~  =  F  (l  -  ygcjps  ig>  m) .  (84) 

Note  that  (83)  can  be  considered  a  version  of  the  Bilby- 
Kroner  decomposition  proposed  for  elastic-plastic  solids  [8, 
62]  and  analyzed  at  length  in  [5, 11,  26]  from  perspectives  of 
differential  geometry  of  anholonomic  space  (neither  A  nor 
S  is  necessarily  integrable  to  a  vector  field).  The  following 
energy  potentials,  measured  per  unit  reference  volume,  are 
defined: 


+  (85) 

W  ^W[A{¥,r])],  /(?7,D)  ^  arf  {l  -  rjf  +  (D)  . 

(86) 

Here,  W  is  the  strain  energy  density  that  depends  on  elastic 
lattice  deformation  A  (assuming  -  =  0  in  (86)), 

/  is  the  total  interfacial  energy  that  includes  a  double-well 
function  with  constant  coefficient  a,  and  =  KabD^D^ 
is  the  square  of  the  fundamental  Finsler  function  given  in 
(76).  For  isotropic  twin  boundary  energy  Kab  =  k^Sab^  and 
equilibrium  surface  energy  Fg  and  regularization  width  Ig 
obey  Kg  =  (3/4)rg/g  and  a  =  12rg//g  [53] .  The  total  energy 
potential  per  unit  volume  is  y/.  Let  O  be  a  compact  domain 
of  Z,  which  can  be  identified  as  a  region  of  the  material  body, 
and  define  the  total  potential  energy  functional 

T[x(X),fj(X)]  =  [  fdn, 

Jn  (87^ 

dfl  =  |det  {Gab)\^^^  dx'  ■  ■  •  a  dX”  a  c/d'  •  •  ■  a  dD^. 

Recall  that  for  solid  bodies,  n  -  p  -  3.  For  quasi-static  con¬ 
ditions  (null  kinetic  energy),  the  Lagrangian  energy  density 
is  2  =  -y/.  The  null  first  variation  of  the  action  integral, 
appropriate  for  essential  (Dirichlet)  boundary  conditions  on 
boundary  dCl,  is 


SI  (Cl)  =  = 


(88) 


where  the  first  variation  of  potential  energy  density  is  defined 
here  by  varying  x  and  rj  within  O,  holding  reference  coordi¬ 
nates  X  and  reference  volume  form  dO  fixed: 


Sy/  - 


dW  dSx‘‘ 


-I- 


/  dW  A  3/ 

\  dr]  dt])  ^  dD^  dX^ 


(89) 


Application  of  the  divergence  theorem,  here  with  vanishing 
variations  of  x  and  rj  on  dCl,  leads  to  the  Euler-Lagrange 
equations  [53] 


<p[,j(X)]  =3[^(X)]'-2[,7(X)]^  (82) 

The  horizontal  part  of  the  deformation  gradient  in  (79)  obeys 
a  multiplicative  decomposition: 

K  =  A’A  (83) 


dx^  dx^ 
df  _  ^  _  dD^  _  ^  _  d^r] 
dr]  ~  dXs  ~  dXAdXs ' 


F  =  A3, 
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The  first  Piola-KirchhofF  stress  tensor  is  P  =  dWjdV,  and 
the  elastic  driving  force  for  twinning  is  c  =  dWIdrj.  These 
equations,  which  specify  mechanical  and  phase  equilibria,  are 
identical  to  those  derived  in  [53] ,  but  have  arrived  here  via 
use  of  Finsler  geometry  on  fiber  bundle  (.  In  order  to  achieve 
such  correspondence,  simplifications  L{X,  D)  L{D)  and 
x{X,D)  —>■  x{X)  have  been  applied,  the  first  reducing 
the  fundamental  Finsler  function  to  one  of  Minkowskian 
geometry  to  describe  energetics  of  twinning  in  an  initially 
homogeneous  single  crystal  body  =  constant). 

A  more  general  and  potentially  powerful  approach  would 
be  to  generalize  fundamental  function  L  and  deformed 
coordinates  x  to  allow  for  all  possible  degrees-of-freedom. 
For  example,  a  fundamental  function  L{X,  D)  corresponding 
to  nonuniform  values  of  Kj^g(X)  in  the  vicinity  of  grain 
or  phase  boundaries,  wherein  properties  change  rapidly 
with  position  X,  could  be  used  instead  of  a  Minkowskian 
(position-independent)  fundamental  function  L(D).  Such 
generalization  would  lead  to  enriched  kinematics  and  nonva¬ 
nishing  connection  coefficients,  and  it  may  yield  new  physical 
and  mathematical  insight  into  equilibrium  equations — for 
example,  when  expressed  in  terms  of  horizontal  and  vertical 
covariant  derivatives  [30] — used  to  describe  mechanics  of 
interfaces  and  heterogeneities  such  as  inclusions  or  other 
defects.  Further  study,  to  be  pursued  in  the  future,  is  needed 
to  relate  such  a  general  geometric  description  to  physical 
processes  in  real  heterogeneous  materials.  An  analogous 
theoretical  description  could  be  derived  straightforwardly  to 
describe  stress-induced  amorphization  or  cleavage  fracture 
in  crystalline  solids,  extending  existing  phase  field  models 
[54,  60]  of  such  phenomena. 

5.  Conclusions 

Finsler  geometry  and  its  prior  applications  towards  contin¬ 
uum  physics  of  materials  with  microstructure  have  been 
reviewed.  A  new  theory,  in  general  considering  a  deformable 
vector  bundle  of  Finsler  character,  has  been  posited,  wherein 
the  director  vector  of  Finsler  space  is  associated  with  a 
gradient  of  a  scalar  order  parameter.  It  has  been  shown 
how  a  particular  version  of  the  new  theory  (Minkowskian 
geometry)  can  reproduce  governing  equations  for  phase 
field  modeling  of  twinning  in  initially  homogeneous  single 
crystals.  A  more  general  approach  allowing  the  fundamental 
function  to  depend  explicitly  on  material  coordinates  has 
been  posited  that  would  offer  enriched  description  of  inter¬ 
facial  mechanics  in  polycrystals  or  materials  with  multiple 
phases. 
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